home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / man / man-part1 / cat4 / flts.4 < prev    next >
Text File  |  1999-09-16  |  2KB  |  133 lines

  1.  
  2.  
  3.  
  4. flts(1)                        Scilab Function                        flts(1)
  5.  
  6.  
  7.  
  8.  
  9.  
  10.  
  11. NAME
  12.   flts - time response (discrete time, sampled system)
  13.  
  14. CALLING SEQUENCE
  15.   [y [,x]]=flts(u,sl [,x0])
  16.   [y]=flts(u,sl [,past])
  17.  
  18. PARAMETERS
  19.  
  20.   u         : matrix (input vector)
  21.  
  22.   sl        : list (linear system syslin)
  23.  
  24.   x0        : vector (initial state ; default value=0)
  25.  
  26.   past      : matrix (of the past ; default value=0)
  27.  
  28.   x,y       : matrices (state and output)
  29.  
  30. DESCRIPTION
  31.  
  32.   State-space form:
  33.  
  34.   sl is     a syslin list containing the matrices of the following linear system
  35.  
  36.   sl=syslin('d',A,B,C,D) (see syslin):
  37.     x[t+1] = A x[t] + B u[t]
  38.     y[t] = C x[t] + D u[t]
  39.   or, more generally, if D is a polynomial matrix (p = degree(D(z))) :
  40.  
  41.     D(z)=D_0 + z D_1 + z^2 D_2 +..+ z^p D_p
  42.     y[t] = C x[t] + D_0 u[t] + D_1 u[t+1] +..+ D_[p] u[t+p]
  43.  
  44.  
  45.  
  46.  
  47.     u=[u0,u1,... un]  (input)
  48.     y=[y0,y1,... yn-p]  (output)
  49.     x=x[n-p+1]  (final state, used as x0 at next call to flts)
  50.  
  51.  
  52.  
  53.  
  54.  
  55.   Transfer form:
  56.  
  57.    y=flts(u,sl[,past]).     Here sl is a linear system in transfer matrix
  58.   representation i.e
  59.  
  60.  
  61.   sl=syslin('d',transfer_matrix) (see syslin).
  62.  
  63.       past = [u     ,...,  u   ]
  64.              [ -nd           -1]
  65.              [y     ,...,  y   ]
  66.              [ -nd           -1]
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  
  73.  
  74.   is the matrix of past values of u and y.
  75.  
  76.   nd is the maximum of degrees of lcm's of each row of the denominator matrix
  77.   of sl.
  78.   u=[u0 u1 ... un]  (input)
  79.   y=[y0 y1 ... yn]  (output)
  80.   p is the difference between maximum degree of numerator and maximum degree
  81.   of denominator
  82.  
  83. EXAMPLE
  84.   sl=syslin('d',1,1,1);u=1:10;
  85.   y=flts(u,sl); [y1,x1]=flts(u(1:5),sl);y2=flts(u(6:10),sl,x1);
  86.   y-[y1,y2]
  87.  
  88.   //With polynomial D:
  89.   z=poly(0,'z');sl=syslin('d',1,1,1,1+z+z^2);p=2;
  90.   y=flts(u,sl);[y1,x1]=flts(u(1:5),sl);
  91.   y2=flts(u(5-p+1:10),sl,x1);  // (update)
  92.   y-[y1,y2]
  93.  
  94.   //Delay (transfer form): flts(u,1/z)
  95.  
  96. SEE ALSO
  97.   ltitr, dsimul, rtitr
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.